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Abstract 

In this work, we apply the anholonomic deformation method for 
constructing new classes of anisotropic cosmological solutions in Ein- 
stein gravity and/or generalizations with nonholonomic variables. The- 
£N1 ' re are analyzed four types of, in general, inhomogeneous metrics, de- 

fined with respect to anholonomic frames and their main geometric 
properties. Such spacetimes contain as particular cases certain confor- 
| mal and/or frame transforms of the well known Friedman-Robertson- 

■ Walker, Bianchi, Kasner and Godel universes and define a great va- 
riety of cosmological models with generic off-diagonal metrics, local 

■ anisotropy and inhomogeneity. It is shown that certain nonholonomic 
| gravitational configurations may mimic de Sitter like inflation scenarios 

and different anisotropic modifications without satisfying any classical 
false-vacuum equation of state. Finally, we speculate on perspectives 
when such off-diagonal solutions can be related to dark energy and 
r> ' dark matter problems in modern cosmology. 
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1 Introduction 



Modifications of general relativity (GR) theory and new classes of cos- 
mological solutions have received much attention as attempts to account 
for dark energy and recent observations from the Wilkinson Microwave 
Anisotropic Probe (WMAP), see [H El 13 H [6] for reviews. There is 
certain evidence of relatively small anisotropic departures from the stan- 
dard Friedmann-Robertson-Walker (FRW) model. However, we argue that 
it may be possible to involve more general classes of anisotropic and/or in- 
homogeneous cosmological solutions, described by generic off-diagonal met- 
ric^] in GR, in order to explain in the bulk the existing experimental data 
and examine various types of cosmological scenaria. 

A series of Bianchi models with anisotropics has been analyzed. There 
is a classification of Bianchi metrics arranging all possible spatially homoge- 
neous models depending on the symmetry properties of their spatial hyper- 
surfaces, see recent developments and references in [3J8]. A few of Bianchi 
universes contain the FRW model as a limiting caseo But even the best 
such models (the so-called Bianchi VII^ class) seem to be inconsistent [9] 
with WMAP data. 

Another classes of anisotropic and/or inhomogeneous solutions are used 
in Kasner and Godel cosmological models [TOj HH [T21 EH HH El HE] . Never- 
theless, the problem to construct the cosmological model and related solu- 
tions of gravitational and matter field equations which would describe most 
realistically the existing experimental data is still unsolved in modern grav- 
ity, cosmology and astrophysics. 

Recently, we developed the so-called anholonomic deformation method 
of constructing exact solutions in gravity |17| (see examples and reviews in 
[HI [181 EQl EH E2] ) • Perhaps, this is the unique existing at present geometric 
method providing a formalism for generating very general classes of solutions 
of gravitational field equations in GR and various high dimension/ metric- 
affine, Lagrange-Finsler, noncommutative or other extensions. The method 
is based on nonlinear connection geometry originally elaborated for Finsler 
spaces, see details in [23j and, in relation to standard theories of physics, in 
[22} 120] . In this work, we shall not concern any questions related to details 
for scenaria of " Finsler cosmology" but consider some geometric techniques 
for generating cosmological solutions parametrized by generic off-diagonal 

1 such metrics can not be diagonalized by any coordinate transforms 

2 The paradigm of modern cosmology is based on the Friedmann-Robertson-Walker 

(FRW) metric which is derived as a spherical symmetric solution of the Einstein equations 

assuming homogeneity and isotropy on large scales. 



3 



metrics. The goal of this paper is to construct and analyze new classes of 
cosmological metrics for four dimensional, 4-d, pseudo-Riemannian space- 
times. 

The off-diagonal anisotropic cosmological solutions to be considered in 
this work are for GR. We shall use an auxiliary linear connection (the so- 
called canonical distinguished connection, d-connection) , for which the Ein- 
stein equations can be solved in general form with respect to some adapted 
nonholonomic frames of reference]^ Such a connection is also metric compat- 
ible and completely defined by the metric structure but contains nontrivial 
torsion terms induced nonholonomically by some off-diagonal coefficients of 
metric. Having constructed generalized exact (cosmological) solutions for an 
auxiliary connection, we have to impose certain additional constraints on co- 
efficients of metrics in order to " extract" exact solutions for the Levi-Civita 
connection. 

We emphasize that the anholonomic deformation method allows us to 
generate "almost all" classes of solutions when the time like coordinate is 
contained as a "nonholonomic" one. Such generic off-diagonal solutions pos- 
sess at least a Killing vector symmetry and, in general, depend additionally 
on two space coordinates. Various models of nonhomogeneous and locally 
anisotropic cosmological models^ can be elaborated. The Einstein equa- 
tions are very complex systems of nonlinear partial differential equations. If 
we perform cosmological approximations for an explicit metric ansatz (for 
instance, considering only the dependence on time by averaging on space 
coordinates and imposing certain additional spacetime symmetries), we get 
some systems of nonlinear ordinary differential equations. Even we may be 
able to solve such systems in a general form, we get only a very restricted 
subclass of cosmological solutions. This way we "cut" the bulk of nonlin- 
ear gravitational interactions and loose a number of important off-diagonal 

3 In modern literature on mathematics and physics, there are considered three equiva- 
lent terms: nonholonomic, anholonomic and/or non-integrable (for convenience, we shall 
use all such terms). Here, we also note that our approach should not be confused with the 
so-called Cartan's moving frame method when some geometric/ physical objects are re- 
defined with respect to some more convenient frames of references/local basis. We consider 
nonholonomic deformations of geometric/ physical objects (for instance, deformations of 
the Levi-Civita connection) by imposing non-integrable constraints on the dynamics of 
gravitational fields and anholonomic frames with associated nonlinear connections struc- 
ture (the last one being defined as a conventional horizontal (h) and vertical (v) spacetime 
splitting) . 

4 we shall use also the equivalent terms 'inhomogeneous / nonhomogeneous' and say 
that the solutions are (locally) anisotropic if the geometric constructions are defined for 
off-diagonal metrics with dependencies both on time and space coordinates 
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solutions. 

In our approach, we can integrate the Einstein equations in some very 
general forms. Performing approximations for general solutions (and not for 
certain coefficients of systems of nonlinear equations) , we find new classes of 
cosmological metrics which extend the already known families of metrics for 
Bianchi universes, Kasner spacetimes etc. We note, that it is not possible to 
derive such cosmological metrics working directly with the Levi-Civita con- 
nection and local coordinate frames. The surprising property of separation 
of equation^ exists for a more general type of connections which are also 
completely (and uniquely, but with a different geometric meaning) deter- 
mined by the metric structure. At the first step, we can construct solutions 
for generalized connections and then, the second step, we have to constrain 
some coefficients of metrics in order to generate (in our case, cosmological) 
solutions in Einstein gravity. Solutions with " un-constrained" off-diagonal 
metrics also present a substantial interest for modern cosmology because 
they can be related to more general models of string/ brane cosmology etc. 

In our further works, we are going to provide an exhaustive study of 
generic off-diagonal and locally anisotropic metrics, and related cosmolog- 
ical models, in Einstein and Lagrange-Finsler theories of gravity that pos- 
sess a FRW limit. The purpose of such constructions is to characterize as 
full as possible the cosmological spacetimes with anisotropics and generic 
off-diagonal nonlinear interactions and thus provide the strongest possible 
constraints on exotic cosmologies. In the present paper, however, we focus 
on the very specific question of whether general off-diagonal cosmological 
solutions can be constructed in Einstein gravity and if such models neces- 
sarily involve de Sitter stadia and possible inflation induced by nonlinear 
gravitational interactions. 

This paper is organized as follows: 

In section [21 we present necessary geometric preliminaries on the nonlin- 
ear connection formalism and nonholonomic deformations of metrics, con- 
nections and frames in (pseudo) Riemannian spacetimes. We outline certain 
classes of important cosmological solutions which in this work will be de- 
formed nonholonomically into generic off-diagonal solutions. The Einstein 
equations are equivalently formulated for two types of important linear con- 
nections. We use the fact that for the so-called canonical distinguished 
connection, the gravitational field equations can be separated with respect 
to adapted nonholonomic frames. 

In section [3l we prove that the Einstein equations can be integrated 

it should be not confused with separation of variables 
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in very general forms containing all possible inhomogeneous and locally 
anisotropic cosmological solutions. There are analyzed some important 
parametrizations and subclasses of such off-diagonal solutions. 

Section U] is devoted to explicit constructions of generic off-diagonal cos- 
mological solutions. We derive families of anisotropic spacetimes containing 
in certain limits the FRW, Bianchi, Kasner and Godel type configurations. 
We show that imposing nonintegrable (nonholonomic) constraints on the 
nonlinear dynamics of off-diagonal gravitational interactions we can model 
various types of anisotropic and de Sitter solutions. 

We summarize and discuss the results in section [5j 

2 Nonholonomic Deformations of Cosmological 
Solutions 

In this section, we give general features of the geometry ofnonholonomic 
deformations and apply this formalism for constructing exact off-diagonal 
cosmological solutions. We follow the notations of [20l [191 [22l [21] were 
details and references can be found. 

2.1 Geometric preliminaries 

Let us consider a (pseudo) Riemannian 4-d manifold V endowed with 
a metric g = g a p{u l )du a ® du" of signature (+,+,—,+) when local co- 
ordinates are parametrized in the form u a = {x l ,y a ), where x % = (x 1 ,^ 2 ) 
and y a = (y 3 = t, y 4 = y) Indices i,j,k,... = 1,2 and a,b, c, ... = 3,4 
are used for a conventional (2 + 2)-splitting of dimension and general ab- 
stract/coordinate indices when a, f3, . . . run values 1, 2, 3, 4. 

We denote by V = {r ^ } the Levi-Civita connection^ with coefficients 
stated with respect to an arbitrary local frame basis e a = (ej,e ) and its 
dual basis e@ = (e J , e h ). Contracting the first and third coefficients of the 
Riemannian curvature tensor 7Z = {R a p-y^} of V, we define the Ricci tensor, 
TZic = {R p$ = R a p a g}, which (in its turn) can be used for computing the 
scalar curvature R = g@ s R where g@ s is inverse to g a p. The Einstein 

6 In our works, we use conventions from [201 122] when left up/low indices are used as 
labels for spaces and geometric objects. We state that y 3 = t because such a parametriza- 
tion will allow us to construct and write down the formulas for equations and solutions in 
a "most" simplified form. 

7 which is uniquely defined by a given tensor g to be metric compatible, Vg = 0, and 
with zero torsion; we follow the conventions established in our previous works [201 [22] , 
including summarizing on "up- low" repeating indices if the contrary is not stated 
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equations on V, for an energy-momentum source of matter T a p, are written 
in the form 

R ps - ^dpsR = xTps, (l) 

where x = const. 

Our goal is to construct exact solutions of gravitational field equations 
([TJ parametrized in the form 

"g = r, t {x\t) °g i {x k : t)dx i ®dx i + r la {x k ,t) °h a {x k ,t)e a ®e a , (2) 
e 3 = dt + V f{x k , t) ° Wl {x k ,t)dx\ e 4 = dy 4 + r]f(x k ,t) ° ni (x k ,t)dx i , 

for certain classes of coefficients (functions) to be defined below. In brief, 
we shall write such metrics as 

g=g tj dx l ® dx j + h ab (dy a + N%dx k )®{dy b + N b k dx k ), (3) 

where, for gij = diag[gi = rn °gi] and h ab = diag[h a = r] a °h a ] 
and iV| = Wi = rjf °Wi and A^ = rii = rjf °rij. The gravitational 'po- 
larizations' ??q. and 77? determine nonholonomic deformations of metrics, 
°g=[ °h a , °N%] — > ,? g=[ gi,h a ,N%}. Such transforms (with deforma- 
tions of the frame, metric, connections and other fundamental geometric 
structures) are more general than those considered for the Cartan's moving 
frame method, when the geometric objects are re-defined equivalently with 
respect to necessary systems of reference. 

Any set of coefficients A^ in ([3]) state on V some A^-adapted frame, 
e a , and dual frame, e' 3 , structures (i.e. N-elongated partial derivatives, 
respectively, differentials) 

e a = Ui = 8i- N?d a , e b = d b = -^j , (4) 

= [j = dx\e a = dy a + N^dx { ) . (5) 
Such local bases satisfy some nonholonomic relations 

T 1 

a/3 

with nontrivial anholonomy coefficients w2 (u) , 



*,ep] = e a ep - e^e^ = w 7 (u) e 7 , (6) 



W V = " W %- = n ij = e 3 N i " <a = ~™ b ai = daNl (7) 

We can fix on V such systems of N-elongated frames when the sets of co- 
efficients N = {N£} define a Whitney splitting (in general, non-integrable) 
of tangent space TV to V, 

TV = hV®vV (8) 
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into conventional horizontal (h) and vertical (v) subspaces, respectively, 
KV and uV. Such a geometric object defines a nonlinear connection (N- 
connection) structure. 

For any metric g J3} on a spacetime V, there is an infinite number of 
metric compatible linear connections D, satisfying the conditions Dg = 0, 
and completely defined by g. A subclass of such linear connections can be 
adapted to a chosen N-connection structure N, when the splitting ([8]) is 
preserved under parallelism, and called distinguished connections (in brief, 
d-connections) . A general d-connection is denoted by a boldface symbol 
D = (hD, vD) distinguished into, respectively, h- and v-covariant deriva- 
tives, hD and vD. To construct exact solutions in gravity theories is con- 
venient to work with the so-called canonical d-connection, D = {r 7 Q ^}, 
which with respect to N-adapted bases (|4|) and © is given by coefficients 
f ^ = (%,L^,q c ,^), for hD = {D jk ,L a bk } and vD = {C^C^}, 
where 

L)k = 2 gir ( e k9jr + ejgkr - e r g jk ) , (9) 
L a bk = e b {N^) + hi ac [e k h bc -h dc e b Ni-h db e c Ni 
C% = \g ik e c 9jk, C£ = \h ad (e c h bd + e c h cd - e d h bc ) . 

This canonical d-connection D and its torsion T = {T 7 a/3 = T 7 a/3 — T 7 ^; 
T l j k ,T l j a ,T a j i ,T a bi ,T° bc }, where the nontrivial coefficients 

T^i j'i j'i i-ni _ p-ii i-pa (\a /'in\ 

jk — ^jk ~ ^kji 1 ja — <^j b i 1 ji ~ ~ i2 jii 

Taj = L c a j — e a [Nj),T a bc = C bc — C" 6 , 

are completely defined by the coefficients of metric g © following the condi- 
tions that Dg =0 and the "pure" horizontal and vertical torsion coefficients 
are zero, i. e. T l - k = and T a bc = 0. 

Any geometric construction for the canonical d-connection D can be 
re-defined equivalently into a similar one with the Levi-Civita connection 
following formula 

where the distortion tensor Z^ a a is constructed in a unique form from the 
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Zbk. — +r 'cd Tkb, %kb — ^jkhcbgi 1 + E % jf, C[ b , Zj k — 0, (12) 



ryil _ — >— 'dd rpC rya r\ ryi 9 ^ 

^jb — ~ ^cb 1 jdi ^bc — u ' ^ab ~ ~~ 



coefficients of a metric g a p, 

711 r<i „ uab 1 r>a yi ^ r\c i _ji i-ji/i pij 

^jk — —^jb9ik n — g ll jki ^bk — 2* L jk n cbg ~ ^jk u hb^ 

2 

2 Tjahcb + Tj b h C a 

for E? k = ±(6*6% -g jk g ih ) and ±S^ = \{5«8 b d + h*h!*). 

2.2 Limits to known cosmological solutions 

In this work, we shall construct new classes of cosmological solutions^ 
with metrics ''g, i.e. 'target' metrics, possessing certain limits, for rj a , 
r]f 1, or ijf —7-0, to °g (a 'prime' metric) which is a conformal, or 
frame/coordinate, transform of a well known FRW, Bianchi, Kasner, or 
another type solution, see reviews of results in [71 El 124"! [25]. 

2.2.1 FRW metrics 

The FRW cosmological solution can be written in the form 



Fg=a 2 (t) ( rj + r 2 d0®d0 ) - dt®dt + a 2 (t)r 2 sin 2 9d<p®dip, (13) 

1 KjT 



with k = ±1,0, when the coordinates and coefficients of metric are paramet- 
rized, respectively, in the form x 1 = r,x 2 = 9, y s = t,y 4 = (p (for spherical 
coordinates) and f9i = o- 2 /(l — nr 2 ), pg2 = a 2 r 2 /(l — nr 2 ), ph^ = 
— 1, i?/i4 = a 2 (t)r 2 sin 2 6 and fN? = |H1 This metric is an exact solution 
of equations ([I]) with a perfect fluid stress-energy tensor, 

T a p = diag[-p, -p, p, -p], (14) 

where p and p are the proper energy density and pressure in the fluid rest 
frame. The Einstein equations for ansatz (I13p take the form of two coupled 
nonlinear ordinary differential equations (also called the Friedmann equa- 
tions) 

9 f a*\ 2 1 k , „. 

V a I 6 a z 



in general, such spacetimes inhomogeneous and anisotropic 
9 Instead of FRW as a 'prime' metric °g, we can consider any Biachi, Kasner etc 
cosmological solutions outlined in [8] [9| HI [5] |6] and references therein. 
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and 

H* + H 2 = — = - l -{p + ?,p)i (16) 
a o 

where the strong energy conditions for matter, p + 3p > 0, must be satisfied 
for an expanding universe^ The Hubble parameter H = ^ has the unit of 
inverse time and is positive (negative) for an expanding (collapsing) universe. 
The equations (|15p and (|16p may be combined (which is also related to the 
condition V a T a ^ = 0) 

p* + 3H(p + p) = 0. 
For simplicity, we can consider k = with 

F g=a 2 (t) (dx(&dx + dz<g>dz) - dt®dt + a 2 (t)dy<g)dy, (17) 

for Cartezian coordinates and coefficients of metric parametrized, respec- 
tively, in the form x 1 = x,x 2 = z, y 3 = t,y 4 = y and °g\ = f92 = 
a 2 , i?/i3 = — 1 and f^a = a 2 (in this case, the nontrivial coefficients of 
metric depend only on time like coordinate, t, but not on space like ones). 

2.2.2 Bianchi type metrics 

All possible spatially homogeneous but anisotropic relativistic cosmolog- 
ical models are arranged following the Bianchi classification corresponding 
to symmetry properties of their spatial hypersurfaces [261 E3 E] • Such cos- 
mologocal metrics can be parametrized by orthonormal tetrad (vierbein) 
bases e a " = e^^/d/du— , when 

Bga"(3" = Bd- a ii Bd-pn B9ap = diag[l, 1, -1, 1] (18) 

and 

i BV, b^"] = BW 1 ^,^, (t) sey, 
with time dependent 'structure constants' 

B w a"/3" (*) = e a"fS"T"n T 7 (t) + 5j„b a ii (t) - P a ,,bpn (t) (19) 

(determined by some diagonal tensor, n T 7 , and vector, b a n, fields) used for 
the mentioned classification. Depending on parametrization of such tensor 

10 For our purposes, and following our former works on geometric methods in gravity 
and exact solutions [171 1191 [20] . it is convenient to introduce a system of notations which 
is different form those used in standard books on cosmology (see, for instance, [241 1251 |4l 
El El 03 [H El i were readers may consult details and references on modern cosmology). 
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and vector objects, there are constructed the so-called Bianchi universes 
which are closed, or not, to the homogeneous and isotropic FRW case. With 
nontrivial limits to observable cosmology, there are the so -called Bianchi 
/, V, VIIq, Vllh and IX cosmologies. 



2.2.3 Kasner type metrics 

For instance, in four dimensional gravity, such a metric is written in the 
form 

K g=t 2pi dx(g)dx + t 2p3 dz(g)dz - dt®dt + t 2p2 dy(g>dy, (20) 

with K gi = t 2pi , K g 2 = t 2p:i , K h 3 = -1, K h 4 = t 2p2 , K Nf = 0, see details 
(including modern brane generalizations) and references in |10t \TT\ [T2] . The 
constants pi,P2,P3 define solutions of the Einstein equations if there are 
satisfied the conditions 

2 3p _ 2p_ ip ; pi) 
L P) ={pi) +(P2) +{P3) , P = P1+P2+P3, P = P1P2+P2P3+P1P3- 

Following the anholonomic deformation method, we shall generalize such 
solutions to generic off-diagonal cosmological configurations, see section I4TT1 



2.2.4 Godel model 

The theoretical model for the study of rotating cosmology is determined 
by the Godel solution [13] 



Gg= go 2 



e 2x 



dx®dx + —^-dz(&dz — (dt — e x dz)®{dt — e x dz) + dy®dy 



(22) 

when G gi = G a 2 , G g 2 = go 2 , = - G« 2 , G^4 = aa 2 , 

G Nf = QWi = —e x , cNf = 0. The matter is described as a dust with the 
energy density g £ an d there is a nontrivial negative (with opposite sign to 
that introduced by Einstein) cosmological constant gA determining the an- 
gular velocity gw of the cosmic rotation, quj 2 = 1/2 qcl 2 = 4irG g£ = 
— gA, with G as Newton's gravitational constant. The parametrizations 
of coordinates and coefficients of metrics are different from that considered 
in the former works, see a comprehencive review of results and references 
in |144 [151 fl6] . In section 14.1} we shall provide an off-diagonal generaliza- 
tion defining rotating universes with polarized cosmological constants and 
nonholonomic/constrained rotations in nontrivial backgrounds. 
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2.3 The Einstein equations for connections D and V 

The Einstein equations (pQ) for a metric can be rewritten equivalently 
using the canonical d-connection D, 

R ps - \m S R = (23) 

L c aj = e a (Nf), q b = 0, n a ji = 0, (24) 

where R p§ is the Ricci tensor for r 7 ^, S R = g^ 5 R ps and T pg is such way 
constructed that Ypg ~ * for D — > V. We emphasize here that if the 

constraints (|2"^j) are satisfied the tensors T 7 ^ ([TUj) and Z 7 ^„ (fT2|) are zero. 
This states that r 7 Qj g = T~ t a p, with respect to N-adapted frames (jU and 

©, see (HU, even D / V. 

In a series of our works [TUJ [2D1 EU [221 E], we provided detailed proofs 
that for constructing exact solutions with generic off-diagonal metrics it is 
more convenient to work with the canonical d-connection D than with the 
Levi-Civita connection V; the last one is considered to be the standard one 
for general relativity. The surprising thing is that the "nonholonomic" grav- 
itational field equations ([23]) split in such a form, with respect to N-adapted 
frames (j4|) and ([5]) , that the resulting system of partial differential equations 
(see below, for instance, (|26p - (|29p ) can be solved in very general forms. 
In order to generate exact solutions for V, we have to impose additional 
constraints (I24D on coefficients of metric g ([3]) (for instance, on ''g ([2]) if 
we wont to generate new classes of cosmological solutions) l n l 

For an ansatz of type ([2]), the Einstein equations (|23h for D with a 
general source of tvpa 12 ! 

T a p = diag[T 7 ; T 1 = T 2 = T 2 (x k ,t); T 3 = T 4 = T 4 (x fc )] (25) 

transform into a system of nonlinear partial differential equations with sep- 
aration of equations for h- and v-components of metric and N-connection 



11 We note that it is not possible to solve the Einstein equations in general form working 
directly with V because this way we do not get a separation of nonlinear partial differential 
equations. Our idea is to use a more general connection (also defined completely by the 
same metric), when the system of nonlinear equations can be integrated in some general 
forms, and than to constrain the solutions to generate metrics for the general relativity 
theory. 

12 such parametrizations of energy-momentum tensors are possible by corresponding 
nonholonomic frame and/or coordinate frame transform for very general matter sources, 
including some important cases with cosmological constants and various models of locally 
anisotropic fluid/scalar field/ spinor/ gauge fields interactions on curved spaces 
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coefficients, 



R2 



25152 



9" 



2gi 2g 2 



I92 2 91 



(26) 



-T 4 (* fe ), 



r\ 



2/13/14 



h\* 



(Kf KK 

2h 4 2h 3 



-r 2 (x k ,t), 



(27) 



R 



Wk 
2/l 4 

4/i 4 
/'1 



+ 



2/i 4 2/i 3 



d k h 



fc«4 



2/t, 



n k + 



7 * 1 * 

h* - -h A 



2/14 

2/i 3 



0, 



0. 



(28) 



(29) 



In brief, we wrote the partial derivatives in the form a* = da/dx 1 , a' = 
da/dx 2 , a* = da/dt. The ansatz ([2]) does not depend on variable y 4 (that 
why we do not have terms with d/dy 4 ). 

The above system of equations can be solved in very general forms for 
arbitrary dimensions and signatures as we proved in Refs. (T9j 1201 |2T | [22 ], [T7] 
(see also next section). In this work, we analyze " cosmological configura- 
tions" for D when y 3 = t for generic off-diagonal metrics of type ([2]). New 
classes of cosmological conditions in general relativity, with the Levi-Civita 
connection V, will be extracted by imposing additional constraints 



w* = ejln \h4\,e k Wi = GiW k , n* = 0, din k = <9 fc n; 
satisfying the conditions 



(30) 



3 General Cosmological Off— Diagonal Solutions 

In this section, we construct in explicit form and analyze the properties 
of possible classes of solutions depending on time like variable t and on 
"horizontal" spacelike coordinates x 1 of gravitational field equations (|26f> — 
(|29p . for the canonical d-connection, and of constraints (|3Q|) selecting Levi- 
Civita configurations. 
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3.1 Type 1: Solutions with h% A ^ and T 2 , 4 ^ 
Such metrics are defined by a metric ansatz 

"g = e^ xk) dx i ^>dx i + h 3 {x k ,t)e 3 ®e 3 + h A {x k ,t)e%e\ (31) 
e 3 = dt + Wi(x k ,t)dx\e A = dy A + ni(x k ,t)dx l 

with the coefficients being solutions of the system^ 3 ! 

4> + iP" = 2T 4 (x fc ), (32) 
h* 4 = 2h 3 h 4 T 2 (x i ,t) /</>*, (33) 

pwi + ai = 0, (34) 

nf+ 7 < = 0, (35) 

where 

= ln| t=M= |, Oi = /»W, P = h\<f>* } j=(ln\h 4 \ 3/2 /\hs\)* . (36) 

For /i| 7^ 0; T 2 ^ 0, we get (ft* 7^ 0. Prescribing any nonconstant eft = 
cj)(x l ,t) as a generating function, we can construct exact solutions of (|32l) 
(|35p solving respectively the two dimensional Laplace equation, for g\ = 
Q2 = e^ 1 ); integrating on t, in order to determine /13, /14 and rii, and 
solving algebraic equations, for Wi. We obtain (computing consequently for 
a chosen (f)(x k ,t)) 

a = » = ^,*. = ± 1*^21, (37) 

1„ 1 2. 



= -di<j)/(t>*, m = 1 n k (x l )+ 2 n k (x l ) J [V(v%J) 

where ®h A (x k ), 1 n k (xM and 2 n/% (x*) are integration functions. We have to 
fix a corresponding sign ± in order to generate a necessary local signature 

of type (+ H h) for some chosen 4>, Y2 and T4. 

Here we note that the general off-diagonal solutions (|13p include as 
particular cases the solutions for a nontrivial cosmological constant Tj = A, 



it is an equivalent of equations (|26|) - (|29[1 for h% 4 7^ 
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or nonholonomic configurations with polarizations of such constants, A — > 
h X(x k ) = T 4 (x fe ) and A v X(x k ,t) = T 2 (x k ,t). 

In order to construct exact solutions for the Levi-Civita connection, we 
have to constrain the coefficients (f3T|) of metric ([3Tj) to satisfy the condi- 
tions (|3"Uj) . Such constraints restrict the class of generating and integration 
functions. For instance, we can chose that 2 n k (x J ) = and 1 n k (x l ) are 
any functions satisfying the conditions d{ 1 rii c = d k 1 nj. The constraints for 
(j)(x k ,t) follow from constraints on N-connection coefficients Wi = —di4>/(j)*, 

(wMT + m[4>] (hMT + 9M<f>] = o, 

di w k [4>\ = d k Wi[4>], (38) 

where, for instance, we denoted by h^cfi] the functional dependence on <fi. 
Such conditions are always satisfied for cosmological solutions with <p = <p{t) 
or if (j) = const (in the last case Wi(x k ,t) can be any functions as follows 
from (fM|) with zero /3 and a^, see (f3Uj) ). 



3.2 Important special cases 

We can construct such solutions for certain special parametrizations of 
coefficients for ansatz (]37p subjected to the condition to be a solution of 
equations (p?2D - (f3"5]) . 



3.2.1 Type 2: Solutions with h\ = 

The equation (p7|) can be solved for such a case, h\ = 0, only if T2 = 0. 
We can consider any functions Wi(x k ,t) as solutions of (j28j) . and its equiv- 
alent (fM|h because the coefficients /3 and «j, see (f36j) . are zero. To find 
nontrivial values of rii we can integrate (135h for h% = and any given 
which results in rii = n k (x J ) + 2 n k [x % \ J h^dt. We can consider any 

9i = 92 = e^ x \ with ^(x k ) determined by (|32p for a given T4(x fe ). 

Summarizing the results, we conclude that this class of solutions is de- 
fined by an ansatz 

= e^W ®dx i + h 3 (x k ,t)e 3 ®e 3 + /i 4 (x fc )e 4 ®e 4 , (39) 



e 3 = dt + Wi(x h ,t)dx\ e 4 = ciy 4 + 



l n k (x l ) + 2 n k (a;*) / fatft 



dx* 



for arbitrary generating functions h^(x k ,t),Wi(x k ,t), °h,4(x k ) and integra- 
tion functions x nfc (x*) and 2 n& (x*) . 
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The conditions (|30p selecting from (|39[) a subclass of solutions for the 
Levi-Civita connection transform into the equations 

2 n k (x 1 ) = and di l nfc = d k l rii, 
w* + di °/i4 = and di w k = d k Wi, 

for any such Wi(x k ,t) and °/i4(x fc ). This class of constraints do not involve 
the generating function h%{x k ,t). 

3.2.2 Type 3: Solutions with h% = and h* 4 ^ 

Such metrics are defined by ansatz of type 

v g = e^ (xk) dx i 0dx i - %(x fc )eW + h 4 (x k ,t)e%e\ 

e 3 = dt + Wi (x k ,t)dx\e 4 = dy 4 + ni{x k ,t)dx\ (40) 

where Q\ = Q2 = \ with ip(x k ) being a solution of (|32p for any given 
T4(x k ). The function h 4 {x k ,t) must satisfy the equation (133 p which for = 
is just 

hT-^-2°h 3 h 4 T 2 (x k ,t) = 0. 
The N— connection coefficients are 

Wi = -dij/f, Hi = V (x l ) + 2 n k {x 1 ) I [l/{yjMf]dt, 



when (f) = In \h* A /sJ\ 

The Levi-Civita configurations for solutions (|40p are selected by the 
conditions (|30p which, for this case, are satisfied if 

2 n fc = and di 1 n k = d k 1 n i , 

and 

(wS\Y + w$] (hS\)* + dihS] = 0, 
di w k [4>] = d k Wi[(j)]. 

Such conditions are similar to (|38p but for a different relation of v-coefficients 
of metric to another type of generating function (j). They are always satisfied 
for cosmological solutions with <j> = <f>{t) or if = const (in the last case 
Wi (x k , t) can be any functions as follows from (|34p with zero /3 and ctj , see 
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3.2.3 Type 4: Solutions with = const 

If we fix (f> = (f>o = const in (fBUj) . but h\ ^ and /i| 7^ 0, we can express 
the general solutions of (|32 p - (|35p in the form 



"g = e^**W ® dz* - °/i 2 [/* (sSt)] 2 |?r (x*,*) 



e 3 (gie 3 



2 ^ e 4 ®e 4 



+/ 2 (**,*) 



e 3 = dt + Wi(x k ,t)dx\ e 4 = dy 4 + n k (x\t) dx i , (41) 

where °h = const, g± = gi = e^ x \ with tp(x k ) being a solution of (|32|) 

for any given T 4 (x fc ), and ?T (x\ t) = ?4 [o] (^) - § / T 2(A *)[/ 2 (z*, *)] 2 <ft. 
The N-connection coefficients Nf = Wi(x k ,t) and Nf = ni(x k ,t) are 

Wi = (42) 

and 

n fc = V (x*) + 2 n fe (**) / ^^41 <* dL ( 43 ) 

We must take qjo] (x l ) = ±1 if sx (x*,i) = ±1 for T2 — > 0. In such a 
case, the functions h 3 = — °h 2 [f* (x*,i)] 2 and /14 = f 2 (x*,t) satisfy the 
equation (|33[) written in the form \/|/i3| = °h(^\h^\)* , which is compatible 
with the condition <p = <pQ. 

The subclass of solutions for the Levi-Civita connection with ansatz of 
type (|4~T1) is subjected additionally to the conditions (|30l) . For instance, we 
can chose that 2 n k (x l ) = and n& (x 4 ) are any functions satisfying the 
conditions di n k = dk 1 ni. The constraints on values Wi = — <9i?x/?x result 
in constraints on ?x, which is determined by T2 and /, 

(u>i[er])* + w i[^r] (M^D* + dih^r] = 0, 

di w k [sr] = d k Wi[<5r], (44) 

where, for instance, we denoted by /ufcr] the functional dependence on qy. 
Such conditions are always satisfied for cosmological solutions with / = f(t). 
For D, if T2 = and 4> = const, the coefficients Wi(x k ,t) can be arbitrary 
functions (we can fix ?x = 1 , which does not impose a functional dependence 
of Wi on ?x) as follows from (j34"l) with zero (3 and on, see (|36j) . To generate 
solutions for V such Wi must be additionally constrained following formulas 
re-written for Wj[?x] —> Wi(x k ,t) and /i 4 [?x] —> ^4 (x*,i) . 
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We note that any solution g = { g a > p> {u a ) } of the Einstein equations (f23|) 
and/or ([I]) with Killing symmetry d/dy (for local coordinates in the form 
y s = t and y 4 = y) can be parametrized in a form derived in this section. 
Using frame transforms of type e a = e a ' a e a >, with g aj g = e a ' a e aga'P'i for any 
Sa/3 ©, we relate the class of such (inhomogeneous) cosmological solutions, 
for instance, to the family of metrics of type (|3ip l 14 l 

Following recent developments from Ref. [T7], we can construct 'non- 
Killing' solutions depending on all coordinates. Such general classes of so- 
lutions can be parametrized in the form 

g = +g i (x k )dx i ®dx i + u 2 (x j ,t,y)h a (x k ,t)e a ®e a , 
e 3 = dy 3 + w i (x h ,t)dx i ,e 4 = dy 4 + ni(x k ,t)dx\ (45) 

for any oj for which e^uj = dkUJ + wt^o* + n^du/dy = 0, when (|45h with 
u) 2 = 1 is of type (2J). 

Finally, we note that the metrics constructed above define general cos- 
mological solutions of Einstein equations for any type of sources xT^ which 
can be parametrized 1 "^! in a formally diagonalized form T 7 (|25p , with respect 
to a nonholonomic frame of reference. 



4 Examples of Cosmological Models with Off— Di- 
agonal Metrics 

The goal of this section is to analyze explicit examples of generic off- 
diagonal cosmological solutions with T4 = but, in general, with non- 
vanishing T2(i). Such solutions are constructed as examples of metrics (|37p. 
(f39j) . (J40j) and ([4"Ijh The new classes of cosmological metrics have respective 
'diagonal' limits to conformal and/or frame transforms of metrics (]17p . (|18p . 
(20}, (J22J. 

We consider a particular parametrization of metrics of type ([2J when 

"g = r, i {x k ,t)°g i (x k ,t)dx i ^dx i + r la (t) h a (t)e a ^e a , (46) 
e 3 = dt + rjf(t) °Wi{t)dx\ e 4 = dy A + r)f(t) ° ni (t)dx\ 

or (for an alternative parametrization) 

e 3 = dt + [r]f(t) + °^(t)] di*, e 4 = dy A + [ry 4 (t) + m(t)] dx\ 

14 We have to solve certain systems of quadratic algebraic equations and define some 
e Q Q,(u' 3 ), choosing a convenient system of coordinates u a = u a (u 13 ). 

15 using chains of nonholonomic frame transforms, this is possible for 'almost' all physi- 
cally important energy-momentum tensors 
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when (for some constructions) g\ = n\(x k ,t) °gi(x k ,t) = 1 and g 2 = 
r]i(x k ,t) °g2{x k ,t) = 1 are trivial solutions of ([26]) . and (f32]h with T4 = 0. 
For homogeneous configurations, we can always introduce such coordinates 
when the coefficients of metrics do not depend on space variables (a formal 
presence of radial and angular variables may exist, for instance, in spherical 
coordinates, like in the case of FRW metric (|13p ) The polarization functions 
ij a (t) and nf(t) can be chosen in a form which allows us to generate homo- 
geneous configurations (i.e. solutions not depending on space coordinates) 
with h a {t) = Va {t) °h a (t) and Wi {t) = nf{t) ° Wi {t), m{t) = r]t(t) ° ni {t). 

4.1 Nonholonomic FRW, Bianchi, Kasner and Godel type 
configurat ions 

4.1.1 N anholonomic FRW generalizations 

Solutions of type 1: We chose °gi = 1, hs(t) = —a~ 2 (t), °/i 4 = 1, 
where a(t) is determined by equations (|15p and (|16p for the FRW model. A 
class of anistoropic and inhomogeneous solutions parametrized by metrics 
of type ([37|) is generated by data 

9i = 92 = 1, n 3 = ± — — ; , 

T 2 (x\t) 

0, ,jfcN , o f(ew[2Hx k ,t)]y 



h 4 = °h A (x k ) ± 2 j 



■dt, 



T 2 (x\t) 

Wi = -di4>/<j)*, m = 1 n k (x i )+ 2 n k (x i ) I [h 3 /(y/\hT\) 3 ]dt. 



Off-diagonal cosmological metrics depending only on variable t can be ex- 
tracted by gravitational polarizations 

tr = l,rj 3 = a\tW(t)\/r 2 (t), m = 1± 2j {eM ^ )]T dt, 

V f(t) = 0, V f(t)= V+ 2 n k [ M*)/( V\mW\f}dt, (47) 



with some x n k = const, 2 n k = const, when Wi(t) = nf(t) + u>i(t) = 0, 
for °Wi{t) = 0, and n^t) = nf(t) + °ni(t), for m(f) = 0. The factor a 2 {t) 
can be included into a generating function (p(t), or into a source T2(i) (we 
can say that it is nonholonomically modelled by such a generating function, 
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alternatively, effective source). The off-diagonal metric 

v g = dx 1 <g> dx l + dx 2 ® dx 2 + ^ ,} dt®dt + (48) 

T 2 {t) 



T a (t) 



[dy 4 + r ?J 4 (t)dx i ]®[dy 4 + ??i 4 (t)dx i ], 



with r]f(t) computed following formulas (|47p defines a class of cosmologi- 
cal solutions of the Einstein equations (|23p for D with a source T"^ = 
diag[Y 7 ; Yi = Y2 = Y 2 (t); Y3 = Y4 = 0]. We can extract Levi-Civita con- 
figurations if we fix, for instance, the integration constants n k = rip. = 0. 
For gravitational polarizations with certain smooth limits 

-a (t), / — — ->• 0, r/i (i) ->• 



T 2 (*) w '7 T 2 (t) 

(this can be satisfied by a corresponding choosing of functions </>(£), Y 2 (i) and 
integration constants rif., 2 nfc), the solutions (|48p transform not just into 
the FRW metric (113p but into its conformal transform (with multiplication 
on factor a~ 2 (t)). Such off-diagonal cosmological solutions have certain sim- 
ilarities (for small nonholonomic deformations) to conformally transformed 
FRW solutions if 0(i),Y 2 (t) are fixed to mimic a(t) as in equations (fT5|) 
and (116p for the Hubble parameter. 



Solutions of type 2: As for the type 1, we also take °<?j = 1, °h^(t) = 
—a~ 2 (t), °/i4 = 1 but generate anistoropic and inhomogeneous solutions 
parametrized by metrics of type (i39j) 



g = dx 1 ® dx 1 + dx 2 dx 2 + /i 3 (x fc , t)e 3 <g>e 3 + °/i 4 (x fc )e 4 ®e 4 , (49) 



„3 j + 1 „,. /k „4 _ 7 ,4 



di + Wi(x K ,t)dx l , e 4 = dy 4 + 



L n fc (z*) + 2 n k {a?) / frgctt 



dx 4 



for arbitrary generating functions h%{x k ,t),Wi{x k ,t), °/i4(x fe ) and integra- 
tion functions 1 n k (x*) and 2 n k (x*) , when Y 2 = 0. The polarization func- 
tions are 

r]i = 1,573 = a 2 {t)h 3 (x k ,t), m = %(x fc ), 

r]f = Wi{x k ,t),7]f = l n k {x i ) + 2 n k (x i ) / h 3 dt, 



with °Wi(t) = and °ni(t) = 0. Such cosmological solutions are anisotropic 
and inhomogeneous and constructed as nonholonomic deformations of a con- 
formally transformed (with multiplication on factor a~ 2 (t)) FRW metric 
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(|13p . We can prescribe polarization functions rj 3 (x k ,t) when /13 = 773 hs(t) — > 
—a~ 2 {t) for 773 — >■ 1. 

The conditions (|30p selecting Levi-Civita configurations transform into 
equations 

2 n fc (V) = and <9j 1 n k = d k 1 n i , (50) 
w* + di hi = and di w k = d k Wi. 

Such constraints can be satisfied for any generating function /13 (x k , t) but 
impose additional constraints on N-coefficients Wi(x k ,t). 

Off-diagonal metrics of type (|49p can be generated to be only with time 
like dependence of coefficients, when /13 = h^{t),Wi = Wi(t) and rtj(t) are 
determined with some constant values of /14, n^, 2 rifc. Such conditions 
are for the Levi-Civita configurations if Wi = const. This defines vacuum 
solutions of the Einstein equations. They transform conformally and non- 
holonomically a FRW universe into certain vacuum Einstein configurations 
which (in this particular case) can be diagonalized by coordinate transforms. 



Solutions of type 3: It is not possible to construct off-diagonal general- 
izations of type (I40p . with coordinate y 3 = t, for FRW universes, when hz = 
— °h 3 (x k ) does not depend on t and °gi = 1, °h 3 (t) = —a~ 2 (t), °h& = 1. For 
such solutions, we can not obtain /13 ~ oT 2 (t) in the limit of trivial polariza- 
tions. Considering inhomogeneous metrics of type (j4"5|) with u = uj(x l ,t,y), 
where a~ 2 (t) can be included into u 2 , we model general inhomogeneous 
solutions obtained via "vertical" conformal transforms and further nonholo- 
nomic deformations. The metrics are of type 

''g = dx 1 <g> dx 1 + dx 2 ® dx 2 + (51) 

a- 2 {t)Zj 2 {x\t,y) \- °h 3 (x k )e 3 ®e 3 + /i 4 (x fc , i)e 4 ®e 4 ] , 

e 3 = dt + w i (x k ,t)dx i ,e 4: = dy 4 + ni(x k ,t)dx\ 

where w = w/a is a solution of d k uj + w k uj* + n k duj/dy = 0. The equation 
(|33p for /ig = transforms into 

hT _ _ 2 o h3hiT2{x k jt) = Q (52) 

Integrating two times on t, we can define hi{x k ,t) for any given °hs(x k ) 
and T 2 (x k ,t). Then we can compute the N-connection coefficients following 
formulas 



-dicP/cj)*, m = l n k (x l ) + 2 n k (a:*) / [1/ ' {yj\h^\f]dt, (53) 
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for <p = \n\h*J^\ °/t 3 /t4||. 

The Levi-Civita configurations are extracted from the set of solutions 
HSU if 

2 n k (x l ) = and di 1 n k = d k 1 n i , 

and 

(w$])* + w$] (hS]f + dih$\ = 0, 

di w k [4>] = d k Wi[(j)}. 

Such conditions are always satisfied for cosmological solutions with = <p(t) 
or if cf> = const (in the last case Wi(x k ,t) can be any functions as follows 
from (pM|) with zero f3 and see in such cases, we must take T2 = 0). 

The first subclass of metrics (|5ip (with solutions depending only on t, 
with nontrivial u(t,y), for constant /i3 and n k ; 2 n k = and nontrivial 
T2(t)) are parametrized in the form 

''g = dx 1 gidx 1 +dx 2 dx 2 +a~ 2 (t)uj 2 (t,y) x (54) 
[- °h 3 dmdt + h A (t)(dy 4 + 1 n i dx i )®(dy 4 + 1 n i dx i )] , 

where /i 4 (t) is any solution of /if - 2 h 3 h A T 2 {t) = 0. 

The second subclass of metrics (|5ip contains arbitrary functions u>j(i) 
but for T 2 (t) = and /14(f) is any solution of /if - (/i|) 2 /(2/i 4 ) = 0. For 
constant °h 3 and 1 n k and 2 rifc = 0, we get 

''g = dx 1 ® dx 1 + dx 2 (g> dx 2 +a~ 2 (t)u} 2 (t,y) x (55) 
[- °h 3 (dt + w i {t)dx i )®{dt + Wi^dx 1 ) 
+hi{t)(dy A + 1 n i dx i )®{dy 4 + V^)]. 

Both types of solutions (154ft and (I55p are for the Levi-Civita connec- 
tion. They can be generalized for the canonical d-connection by introduc- 
ing nontrivial constants 2 n k in the formulas for N-connection coefficients, 
Nf = rii(t), see ([53]) . with /i 4 (i) respectively computed as in equation ([52]) . 
when T2 = ^2(t) and or T2 = 0. 

For the metric (|51|) . the polarization functions are °h 3 (x k ) 

rji = l,n 3 =u 2 (x\t,y)°h 3 (x k ),r ]4 = a- 2 (t)u 2 ( y x\t,y) h 4 (x k ,t), 
ri = Wi {x k ,t),r,t = 1 n k (x i )+ 2 n k (x i ) J h 3 dt, 

with °Wi(t) = and °rii(t) = 0. We have to eliminate respectively the 
dependencies on x % in these formulas in order to compute the gravitational 
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polarizations for the metrics (|54p and (|55p . For simplicity, we omit the 
formulas for this class of solutions. Such classes of metrics limit to the 
conformally transformed FRW one (with conformal factor a~ 2 (t)) if r\ a — > 1 
and 7]f — » 0. 



Solutions of type 4: If we fix 4> = </>o = const in (|36p . but 7^ and 

/14 ^ 0, we have to reparametrize the sets of generating and integration 
functions in a different form. The solutions (|41[) with °gi = 1, °h^(t) = 
—a~ 2 (t), °/i4 = 1 are of type 

T 'g = dx 1 <g) dx 1 + cfx 2 <g> dx 2 

- °h 2 [f* (x\t)] 2 \ n {x\t) |eW + f (x\t) e%e\ 
e 3 = dt + Wi(x k ,t)dx\ e 4 = dy 4 + n k (x\t) dx\ (56) 

where 



ffr = l-||r 2 (x fc ,t)[/ 2 (xV)] 2 <ii, 



n fc = x n fc (*<) + 2 n k (**) / ^ ]' (x\t) dt. 

The Levi-Civita configurations are chosen when 2 n k = and l n k (x l ) 
are any functions with di 1 n k = d k l rii. The constraints for Wi = — dj^rAr 
are similar to (|44p . 

We obtain from (|56p generic off-diagonal homogeneous cosmological so- 
lutions if we consider metrics with dependence only on t. They can be written 
in the form 

= dx 1 ® dx 1 + dx 2 ® dx 2 - °h 2 [f* (t)} 2 \q r (t) |e 3 ®e 3 + f (t) e 4 ®e 4 , 
e 3 = dt + Wi{t)dx\ e A = dy A + n k {t)dx\ (57) 

where, for constant 1 n k , 2 n k and and cr (*) = 1 — jf / T2(i)[/ 2 (i)] 2 ^) 

«>< = -%T (t) At (t) , n k = x n k + 2 n k [ ^~T^l (*) 

For simplicity, we omit the explicit formulas for gravitational polarizations 
(they can be computed similarly to those given for above examples). 
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Taking 2 rik = for ([57|) and imposing constraints of type (}4"4"]) for Wi(t), 
we generate new classes of cosmological solutions in general relativity. For 
trivial gravitational polarizations and vanishing N-connection coefficients, 
such metrics result in the conformal transform (with factor a~ 2 (t)) of the 
FRW metric (TT3T) . 

4.1.2 Locally anisotropic Bianchi spacetimes 

We speculate how a Bianchi metric Bg a "/3" {32]) can be generalized into 
locally anisotropic solutions. Taking any set of coefficients B9ap(t), we have 

to construct certain frame transform, ° B g a i3 = B&~ a n Be^„ r q^r F^l when 
B g a /3 is parametrized as a prime metric 

° B g= °g l dx i ® dx l + °h a {t)e a ®e a , e 3 = dt+ ° Wi {t)dx\ e 4 = dy A + °ni(t)dx\ 

Introducing corresponding gravitational 77-polarizations, we construct non- 
holonomic deformations ^g — > B g, where the target metric can be any one 
of type or (for uj 2 = 1) of type ([2]). In general, such configurations 

are locally anisotropic and inhomogeneous when the solutions depend on all 
coordinates. 

We can derive generic off-diagonal solutions with the coefficients depend- 
ing only on t, or being some constants, as we proved in details in previous 
sections. The main difference is that for trivial gravitational polarizations 
such metrics describe frame/conformal transforms of Bianchi spacetimes and 
not of conformal transforms of FRW metrics. Explicit constructions depend 
on "structure constants" (|19p for the prime metric and are, in general, dif- 
ferent for four types of off-diagonal generalizations. The length of this paper 
does not allow us to provide such details and analyze possible applications 
in cosmology. 

4.1.3 Off diagonal Kasner type metrics 

The data for the primary metric are taken in the form °g\ = 1, °gi = 

i 2(p3-pi) j 0^3 _ -t-Spi^ °/j 4 _ t 2(p 2 -pi) and o N a = q with constants pu p 2 

and P3 as in the metric (I20p . For simplicity, we analyze in this section 
solutions with p% = pi and consider an example when a Kasner universe is 
generalized to locally anisotropic configurations of type 4 characterized by 

16 we have to solve certain quadratic algebraic equations in order to define frame coeffi- 
cients depending on coordinate t, or being some constants 
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gravitational polarizations 

th = 1, V 3= °h 2 [r(x\t)]\ m = f{x\t), 
r)f = Wi {x\t),rf = ni (x k ,t). 

For h a = f] a °h a and Nf = rjf + °iVf , the target metric is of type (|4"Tj) 
generated for T2 = 0, 

v g = dx 1 <g> dx 1 + dx 2 <g> cte 2 (58) 

- °h 2 [f* (x\ t)} 2 r^eW + / 2 (xV) i^e W, 
e 3 = dt + Wi(x h ,t)dx\ e 4 = dy 4 + n k (x\t) dx\ 

where we can consider arbitrary w% = Wi (x l , t) and 

n k = l n k (x l ) + 2 n k {x*) [ ^'ff ' dt. 

J [f(x l ,t)\ 

The coefficients h a are solutions of the Einstein equations for the canonical 
d-connection, see (|2"7j) (and/or, equivalently, the equation written in 

the form y^ft^l = Q h(\f\h^\)* for arbitrary generating function / if 

= Pi- We have to impose additional constraints on / (x\i) if the last 
condition is not satisfied. 

In the limit of trivial polarizations, the metric (|58p results into a confor- 
mal transform (with factor t 2pi ) of Kasner solution (|20p . In general, such a 
prime metric is not a solution of the Einstein equations for the Levi-Civita 
connection but it is possible to chose gravitational polarizations generating 
vacuum off-diagonal Einstein fields even the conditions of type (|2ip are not 
satisfied. Such target metrics may be stable (it is necessary an additional 
analysis of stability properties for any explicit case). 

We can eliminate dependence on space coordinates and generate solu- 
tions of type 

n g = dx 1 <g> dx 1 + dx 2 <g> dx 2 

- °h 2 [f* (t)} 2 t- 2pi e 3 ®e 3 + f (t) t- 2pi e 4 ®e 4 , 
e 3 = dt + Wi(t)dx\ e 4 = dy 4 + n k (t)dx k , 

for arbitrary Wi = Wi(t) and constant l n k and 2 n k , when n k = 1 n k + 
2 n k J dt [f* (t)] 2 I [f (t)] 2 . To extract Levi-Civita configurations we must fix 
2 n k = and impose constraints of type (|44| on Wi(t). 

In a similar form, we can construct nonholonomic deformations of the 
Kasner universes of types 1-3 and/or to generalize them to solutions of type 
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(|45p . The corresponding target metrics may be with "gravitational chaos" 
or constrained nonholonomically to became stable. We omit such details in 
this work. 

4.1.4 Rotating cosmologies with local anisotropy 

We can chose the prime's metric data to be given by the Godel solution 
fl22) when ° gi = G g t {x), °h a = G h a and °N* = G jV«, with local 
coordinates x % = (x,z) and y a = (t,y). Considering polarizations 

m = e^ z \ m = 2e^- 2x , m = V3 (x,z,t), m = m (x,z,t), 
rjf = Wi(x,z,t),7]f = n k (x,z,t) , 

for g { = rji °gi, °h a = r] a °h a , Nf = rjf + °iVf, when Nf = Wi(x,z,t) and 
Nf = rii(x, z, t), we generate metrics of type 

= G a 2 [e* x '^ (dx®dx + dz(g)dz) 

—773(2;, z, t)(dt + Wi(x, z, t)dx l )<g>(dt + Wi(x, z, ^dx 1 ) 
+r]4(x, z, t)(dy + rii(x, z, t)dx l )®{dy + rij(x, z, t)dx 1 )]. (59) 

Choosing a source determined by cosmological constant T2 = T4 = G X, we 
construct a class of type 1 exact solutions if 

$ + v" = 2 G x, nt = 2 G a 2 vm G X/ct>*, 
fiwi + on = 0, n** + 777,* = 0, 

with the coefficients 

77* 

a>i = G a 2 T)ldi<j), f3 = G a 2 rf 4 <f>* , for <f>(x,z,t) = In |— =====|, 

Vl G a 

7 = (ln \ G a 2 i] 4 \ 3/2 /\ G a 2 rj 3 \j . 

Such solutions are derived directly from prime Godel metrics and have limits 
to rotating universes for trivial polarizations. 

In explicit form, we can model locally anisotropic and inhomogeneous 
models with rotation when h\ / 0; for G X 7^ 0, we get (ft* 7^ 0. We obtain 
(computing consequently for a prescribed generating function cft(x,z,t)) 

Vi = e^ z \ m = 2e^- 2x , V3 = ± gX- 1 x \<f>*(x,z,t)\, 

m = °rfi(x,z) ± 2 G A _1 x exp[2 <ft(x, z,i)]dt, (60) 

Wi = -di<f>/4>*, 

m = 1 n k (x,z)+ 2 n k (x,z) [rj 3 (x, z , t) / \r] 4 (x , z,t)\f]dt. 
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The metric (|59p with coefficients (|60p is an explicit example of solutions 
of type (|37|) when the source is determined by a cosmological constant. It 
defines a rotating cosmology additionally imbedded into nontrivial gravita- 
tional backgrounds, which (in general) are locally anisotropic and inhomo- 
geneous. We can impose restrictions of type (PU|) and select Levi-Civita 
configurations. It is possible also to construct models with anisotropic ro- 
tation when such solutions do not depend on x l , but only on t, or with 
generalizations of the Godel model determined by off-diagonal solutions of 
type 2-4. 

4.2 Modeling anisotropic de Sitter configurations 

We can generate cosmological solutions when the coefficients gi(x k ,t) 
in ([2]) depend explicitly on variable t. Let us consider a conformal factor 
q(x k ,t) when 

9 g = q 2 (x k ,t)[r ]i (x k ) °g i dx i ® dx { + 

r) 3 (x k ,t) °/i 3 (i)e 3 ®e 3 + °/i 4 (x fc )e 4 ®e 4 ], 
e 3 = dt + Wi(x k ,t)dx\e 4 = dy 4 + n i (x k ,t)dx i . (61) 

By straightforward computations, we can prove that the Riemann and Ricci 
tensor for an arbitrary metric compatible d-connection, see details in [21] . 
do not change under transform v g= [gij,h a b, N" 1 ] — > q g= [q 2 gij,q 2 h a b, N" 1 ] 
if 

e lQ = dq/dx i - Wl q* = 0. (62) 

For an ansatz of form (I6ip . the Einstein equations for the canonical d- 
connection (j23|) with source (|25p is equivalent to equations (|26p - (|29p and 
additional equations ([62]) for q(x k ,t). 

We search a subclass of inhomogeneous solutions when q = q{x k ,t)a 2 (t), 
with a prime metric °gi = 1, °h 3 = a~ 2 (t), °/i 4 = °h A (x k ), °Nf = 0, and 
77-polarizations chosen such a way that data 

gi = m{x k ), h 3 = m(x k ,t) °h 3 {t), h 4 = °h 4 (x k ), 
Wi = r,!(x k ,t)+ °iV 3 , n l = r ] ?(x k ,t)+ °N 4 , 

generate a metric (solution of the Einstein equations) of type ()49|) with 
arbitrary generating functions Wi and h 3 . For trivial polarizations, when 
r] a — > 1 and rjf — >■ 0, and q = °h 4 = 1, the metric (|6ip transform into 
the FRW metric (|17p if we take a 2 (t) to be determined as a solution of the 
Friedmann equations (fT5l) and (fTBT) . 
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The goal of this section is to prove that by corresponding nonholonomic 
distributions (constraints) on inhomogeneous metrics we can model de Sitter 
like (exponential on t) cosmological solutions of vacuum Einstein equations. 
In such a model, a function a(t) = aoe Ht contains the Hubble constant H as 
an experimental parameter determining a class of nonholonomic constraints 
and related conformal transform with factor q = q(x k , t)aoe Ht , when the 
equations ([621) are 

dln\q\/dx i - 2 Wi H = 0. 

Parametrizing Wi = t Wi(t) + s Wi(x k ) when d k s Wi = d{ s wj c , we solve this 
equation and satisfy the Levi-Civita conditions ([50]) if 

r]f = 1 n i (x k ), with d k 1 n i = c\ 1 n k , 
s Wk = -d k °/i 4 and \n\q{x k ,t)\ = -2H J s w i (t)dx i - <7t 4 (<). 

Putting together the above formulas in (j6ip . we get a class of inhomo- 
geneous off-diagonal cosmological solutions with local anisotropy, 

9g = ale 2Ht exp[-2H J t w i {t)dx i - °h±(x i )] [dx <g> dx + dz ® dz 

+rj 3 (x k ,t)aQ 2 e- 2Ht e 3 ®e 3 + °h A (x k )e A ®e 4 }, (63) 
e 3 = dt+[ t w i (t)+ s Wi(x k )]dx\e A = dy A + 1 n i {x k )dx i . 

This metric models a de Sitter like expansion with arbitrary generating/ 
integration functions ( t Wi(t), s Wi(x k ), r]s{x k , t), °h±{x l ) l rii(x k )) and con- 
stants ao and H which must be chosen following certain bounary /symmetry 
and other physical superpositions to satisfy the experimental data. We 
can fix such polarization functions (i.e. generating/integration functions) 
when a very short accelerated "inflationary" stage is dominate by a locally 
anisotropic and inhomogeneous vacuum solution of the Einstein equations, 
lasting ~ 10~ 36 and containing the de Sitter metric. Then we can say that 
such an locally anisotropic stage was followed by a decelerated homogeniz- 
ing expansion, first with a radiation dominated era and then by matter 
dominated era. 

The class of off-diagonal inhomogeneous solutions (I63p is different for 
the "diagonal" family of Szeres-Szafron metrics considered in inhomoge- 
neous cosmology. Here we note that the books [28l [29] , including references 
within, provide a comprehensive review of the characteristics, properties 
and exact and/or inhomogeneous cosmological solutions. From any such 
diagonal and off-diagonal inhomogeneous metric, we can recover the FRW 
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model, consider solitonic perturbations and analyze contributions of a non- 
trivial cosmological constant }30j. 

Our main conclusion is that using generic off-diagonal exact solutions of 
the Einstein equations, with correspondingly prescribed nonholonomic dis- 
tributions, we can elaborate cosmological models with exponential expan- 
sion and limits to FRW configurations without additional scalar fields which 
would be responsible for inflation. In our approach, inflation is modelled by 
nonlinear off-diagonal interactions and constraints on such a "pure" gravi- 
tational dynamics (in {31 j , we studied a model of anisotropic brane inflation 
with off-diagonal metrics). 

5 Outlook, Discussion and Conclusions 

In this work we have provided the essential features on applications in 
cosmology of the anholonomic deformation method which relate the geom- 
etry of nonholonomic manifolds to exact solutions in gravity. In particular, 
given any physically important cosmological solution of Einstein equations 
(defining respectively the FRW, Bianchi, Kasner, Godel or other universes), 
we constructed new classes of generic off-diagonal exact cosmological so- 
lutions which are, in general, locally anisotropic and/or inhomogeneous. 
Alternately, the nonholonomic deformations described here for cosmological 
solutions can be viewed as examples of the geometry of nonholonomic distri- 
butions and generalized transforms of geometric structures for classical and 
quantum (non) commutative spacetimes. While we have considered the ap- 
proach for Einstein's gravity, it is clear that the general constructions can be 
extended to higher dimensions, various Lagrange-Finsler and string/brane 
gravity models (albeit with increasing computational complexity in comput- 
ing off-diagonal higher dimension terms and/or higher order nonholonomic 
constraints) . 

Our approach allows us to construct general cosmological solutions in 
various gravity theories with metrics g a /3(u T ) parametrized in the form 

<?1 + LU 2 (wi 2 h3 + OJ 2 (n! 2 /l4) CJ 2 (w 1 W 2 /l3 + ni ri2/l4) LU 2 W 1 ll3 lo 2 n-JiA 

2 uj 2 (w 1 w 2 h3 + n 1 n 2 h4) g 2 + uj 2 (w 2 2 hz + n 2 2 114) lo 2 w 2 hs ui 2 n 2 /i4 

uj 2 w 1 h,3 lo 2 w 2 hs /13 

uj 2 n-JiA uj 2 n 2 hn hi 

where the local coordinates are of type u T = (x*,y a ), for x % = (x 3 ,x 2 ) and 
y a = (y 3 = t,y 4 = y) and spacetime signature (+, +, — , +). The coefficients 
gf.(x i ),ha(x l ,t),Wk(x i ,t),rik(x' l ,t),q(x i ,t) and u(x l ,t,y) can be defined in 
explicit form (following well defined and quite simple procedures) by in- 
tegrating and/or differentiating some generating functions. Such metrics 
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depend on certain classes of integration functions and constants in order to 
define very general classes of exact solutions in Einstein gravity and gen- 
eralizations. As a matter of principle, any solution of gravitational field 
equations with certain general matter fields sources can be represented in 
the above generic off-diagonal form by corresponding frame and coordinate 
transforms [TTj. We have to involve certain additional physical considera- 
tions, suppositions on symmetry of interactions and boundary conditions in 
order to model certain realistic cosmological models and scenaria. 

An important issue which we have briefly discussed in this work concerns 
the most general classes of cosmological solutions with " nonohlonomic" time 
like variable y 3 = t. The priority of the anholonomic deformation method is 
that we can determine almost all possible types of cosmological metrics and 
deformations of connection in general form not making approximations with 
any terms in the associated systems of partial differential equations. The 
surprising propriety of the introduces nonholonomic deformations of fun- 
damental geometric and physical objects is that the constructions are such 
way performed that we get separations of the gravitational field equations 
(with respect to certain adapted frames of reference) which allows us to 
generate exact solutions. Further approximations (for instance, for generic 
off-diagonal metrics depending only on time variables t and with certain 
prescribed spacetime symmetries) are possible, but in such case there are 
not lost important types of nonlinear interactions/evolutions which can be 
"lost" during approximations for deriving effective (simplified) systems of 
equations. 

Furthermore, we have considered various generalizations of the known 
important cosmological solutions which may lead to interesting new insights 
into modern cosmology with locally anisotropic and inhomogeneous physi- 
cal scenaria and nonlinear interactions. These provide an interesting arena 
for further exploration in gravity theories and cosmology. In principle it 
is possible to model various nonstandard inflation, dark energy and dark 
matter effects not introducing additional/exotic scalar and other fields but 
imposing certain nonholonomic constraints on the off-diagonal dynamics of 
gravitational interactions in standard Einstein gravity. In addition, it is 
clear that there are many interesting directions that can be studied within 
the framework of the geometry of nonholonomic distributions/frames and 
off-diagonal metrics in gravity and cosmology. Such investigations became 
possible after a general geometric method of constructing exact solutions 
in gravity was elaborated. 
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